Multi-omics characterization of the monkeypox virus infection

Multiple omics analyzes of Vaccinia virus (VACV) infection have defined molecular characteristics of poxvirus biology. However, little is known about the monkeypox (mpox) virus (MPXV) in humans, which has a different disease manifestation despite its high sequence similarity to VACV. Here, we perform an in-depth multi-omics analysis of the transcriptome, proteome, and phosphoproteome signatures of MPXV-infected primary human fibroblasts to gain insights into the virus-host interplay. In addition to expected perturbations of immune-related pathways, we uncover regulation of the HIPPO and TGF-β pathways. We identify dynamic phosphorylation of both host and viral proteins, which suggests that MAPKs are key regulators of differential phosphorylation in MPXV-infected cells. Among the viral proteins, we find dynamic phosphorylation of H5 that influenced the binding of H5 to dsDNA. Our extensive dataset highlights signaling events and hotspots perturbed by MPXV, extending the current knowledge on poxviruses. We use integrated pathway analysis and drug-target prediction approaches to identify potential drug targets that affect virus growth. Functionally, we exemplify the utility of this approach by identifying inhibitors of MTOR, CHUK/IKBKB, and splicing factor kinases with potent antiviral efficacy against MPXV and VACV.

For all statistical analyses, confirm that the following items are present in in the figure legend, table legend, main text, or or Methods section.

n/a Confirmed
The exact sample size (n) for each experimental group/condition, given as as a discrete number and unit of of measurement A statement on on whether measurements were taken from distinct samples or or whether the same sample was measured repeatedly The statistical test(s) used AND whether they are one-or or two-sided Only common tests should be described solely by name; describe more complex techniques in the Methods section.

A description of of all covariates tested
A description of of any assumptions or or corrections, such as as tests of of normality and adjustment for multiple comparisons A full description of of the statistical parameters including central tendency (e.g.means) or or other basic estimates (e.g.regression coefficient) AND variation (e.g. standard deviation) or or associated estimates of of uncertainty (e.g.confidence intervals) For null hypothesis testing, the test statistic (e.g.F, t, r) with confidence intervals, effect sizes, degrees of of freedom and P value noted Give P values as exact values whenever suitable.
For Bayesian analysis, information on on the choice of of priors and Markov chain Monte Carlo settings For hierarchical and complex designs, identification of of the appropriate level for tests and full reporting of of outcomes Estimates of of effect sizes (e.g.Cohen's d, Pearson's r), ), indicating how they were calculated Our web collection on statistics for biologists contains articles on many of the points above.The MaxQuant output files were imported into R with the in-house msimportr R package (https://doi.org/10.5281/zenodo.7746897)and analyzed with the msglm package (https://doi.org/10.5281/zenodo.7752068)apart from one modification: the instrumental error of of MS MS intensities was modeled with a mixture of of Gaussian and Cauchy instead of of Laplacian distribution.The MS1 intensities of of protein group/PTMspecific LC LC peaks (evidence.txttable of of MaxQuant output) were used for modeling the protein group or or PTM intensities.Transcriptomic data normalization, differential expression analysis and p-value adjustment were performed by by the DESeq2 package (version 1.34.0).For gene set enrichment analysis, we we used the in-house Julia package OptEnrichedSetCover.jl (https://doi.org/10.5281/zenodo.4536596).For upstream promoter analysis regulated transcripts were submitted to to ChEA3 web-based application (https://maayanlab.cloud/chea3/).We We performed uniform manifold approximation and projection (UMAP) dimensionality reduction in in R (4.0.2) using R package UMAP104 (0.2.6.0) and manually annotated thus obtained clusters.The prediction of of host kinase motifs on on the host and viral proteins was performed by by using the Kinase Library toolbox (kinaselibrary.phosphosite.org).In In silico prediction of of protein structure was performed using Alphafold (version 2.3.1) and visualized with Chimera (version 1.3).Electrostatic surface potential of of the modeled structure of of MPXV H5 H5 dimer was calculated by by using the PyMOL plugin APBS electrostatics.Molecular nature portfolio | reporting summary

April 2023
For manuscripts utilizing custom algorithms or software that are central to the research but not yet described in published literature, software must be made available to editors and reviewers.We strongly encourage code deposition in a community repository (e.g.GitHub).See the Nature Portfolio guidelines for submitting code & software for further information.

Data Policy information about availability of data
All manuscripts must include a data availability statement.This statement should provide the following information, where applicable: -Accession codes, unique identifiers, or web links for publicly available datasets -A description of any restrictions on data availability -For clinical datasets or third party data, please ensure that the statement adheres to our policy Research involving human participants, their data, or biological material Policy information about studies with human participants or human data.See also policy information about sex, gender (identity/presentation), and sexual orientation and race, ethnicity and racism.

Reporting on sex and gender
Reporting on race, ethnicity, or other socially relevant groupings

Recruitment
Ethics oversight Note that full information on the approval of the study protocol must also be provided in the manuscript.

Field-specific reporting
Please select the one below that is the best fit for your research.If you are not sure, read the appropriate sections before making your selection.

Life sciences Behavioural & social sciences Ecological, evolutionary & environmental sciences
For a reference copy of the document with all sections, see nature.com/documents/nr-reporting-summary-flat.pdfgraphics depictions were produced with the PyMOL software.
In-house R and Julia packages and scripts used for the bioinformatics analysis of the data have been deposited to public GitHub repositories: msglm: https://doi.Please specify the socially constructed or socially relevant categorization variable(s) used in your manuscript and explain why they were used.Please note that such variables should not be used as proxies for other socially constructed/relevant variables (for example, race or ethnicity should not be used as a proxy for socioeconomic status).Provide clear definitions of the relevant terms used, how they were provided (by the participants/respondents, the researchers, or third parties), and the method(s) used to classify people into the different categories (e.g.self-report, census or administrative data, social media data, etc.) Please provide details about how you controlled for confounding variables in your analyses.
Describe the covariate-relevant population characteristics of the human research participants (e.g.age, genotypic information, past and current diagnosis and treatment categories).If you filled out the behavioural & social sciences study design questions and have nothing to add here, write "See above." Describe how participants were recruited.Outline any potential self-selection bias or other biases that may be present and how these are likely to impact results.
Identify the organization(s) that approved the study protocol.

nature portfolio | reporting summary
April 2023

Life sciences study design
All studies must disclose on these points even when the disclosure is negative.The sample sizes were chosen from past knowledge on the good sample size to ensure adequate power.Sample sizes are always indicated in figure legends or related "Methods" section.
No data were excluded For Mass spectrometry, in vitro viral replication experiments, WB analysis, reporter assay and other bioassays, a minimum of three biological experiments were performed independently.All replications were successful.
N/A.No randomization was used given the small number of samples and the lack of influence of randomization on the experimental design and experimental approach used.(no animal experiments were performed in this study).
N/A. Investigators were not blinded to experimental groups (in vitro experiments required prior knowledge for data interpretation).
Antibodies were validated by infection experiments in this study.Additionally, antibodies have been validated by the manufacturers.
The identity of all the immortalized cell lines used in this study was confirmed by STR-profiling (Eurofins Medigenomix).Official certification can be provided upon request.
All cell lines were tested to be mycoplasma free by standard PCR-based assay.
No commonly misidentified cell line was used in this study.
the terms sex (biological attribute) and gender (shaped by social and cultural circumstances) carefully in order to avoid confusing both terms.Indicate if findings apply to only one sex or gender; describe whether sex and gender were considered in study design; whether sex and/or gender was determined based on self-reporting or assigned and methods used.Provide in the source data disaggregated sex and gender data, where this information has been collected, and if consent has been obtained for sharing of individual-level data; provide overall numbers in this Reporting Summary.Please state if this information has not been collected.Report sex-and gender-based analyses where performed, justify reasons for lack of sex-and gender-based analysis.
We require information from authors about some types of materials, experimental systems and methods used in many studies.Here, indicate whether each material, system or method listed is relevant to your study.If you are not sure if a list item applies to your research, read the appropriate section before selecting a response.